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We prove that the upper electronic bands of K-phase BEDT-TTF salts are adequately modeled 
by an half-filled tight-binding lattice with one site per cell. The band parameters are derived from 
recent ab initio calculations, getting a very simple but extremely accurate one-electron picture. This 
picture allows us to solve the BCS gap equation adopting a real-space pairing potential. Comparison 
of the calculated superconducting properties with the experimental data points to isotropic so— 
pairing. Residual many-body or phonon-mediated interactions offer a plausible explanation of the 
d ■ large variety of physical properties observed in K-phase BEDT-TTF salts. 
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Oh' Superconductivity (SC) in organic charge transfer 
(CT) salts has been discovered more than 15 years ago. 
At present, K-phase BEDT-TTF (ET) salts are the most 
C$ ', promising organic superconductors (OSC) OSC have 
similar features to cuprate superconductors: they exhibit 
I , highly bidimensional structures, narrow electronic bands, 
' O ■ low carrier densities and possibly strong electronic cor- 
relations. Superconducting properties are also similar in 
q ■ organics and cuprates, exhibiting singlet pairing very 
low coherence lengths, high magnetic penetration depths 
| ■ and critical fields Q. Perhaps the most characteristic 
feature of cuprate superconductors is the competition be- 

■ tween SC and antiferromagnetism (AFM). This competi- 

\ tion shows up also in OSC, as demonstrated by the pres- 
ence of antiferromagnetic spin fluctuations in the metallic 
state of several compounds ||, or even by reentrant SC 

r-.^, ■ in the presence of magnetic order. In this respect the 
different behavior of k-ET 2 Cu[N(CN) 2 ]X (ET-X) salts 

■ where X=C1, Br is noteworthy: the two compounds are 
isostructural, but at ambient pressure the Br-compound 

Q ■ is a superconductor with T c = 11.6 K, whereas the Cl- 
1 , \ compound is a Mott insulator with (possibly commensu- 
rate) antiferromagnetic ordering ||. On the other hand, 
Q \ even a very small external pressure (270 bar), turns ET- 
O ■ CI to the superconducting ground state (T c = 12.8 K) 
jij . Experimental studies of normal and superconducting 
V "~j ■ state properties in OSC are still scanty and often incon- 
rS | elusive, due to the extreme sensitivity of the material to 

■ pressure, radiation and sample preparation. At the same 
time, the amount of theoretical work on OSC is scarce if 
compared with that developed for cuprates, possibly due 
to the apparent complexity of the OSC structures. In the 
following we will show that, by a proper choice of the ba- 
sis functions, the electronic structure of K-phase ET salts 
can be accurately described in terms of a rectangular lat- 
tice with nearest and next-nearest neighbor interactions. 
Similar one-band models are often adopted in the study 
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of cuprates, however these models lead to an oversim- 
plified view of the complex cuprate band-structures ||. 
OSC then constitute an unique opportunity to test SC 
theories within a simple model, yet accurate enough to 
allow a quantitative comparison with experimental data. 

The basic structural unit of K-phases is a pair of 
ET molecules. These pairs (dimers) arrange themselves 
nearly perpendicularly in a check-board pattern to con- 
struct two dimensional cation layers intercalated by anion 
sheets. As usual for organic CT salts, the relevant physics 
of ET compounds is described in terms of electrons hop- 
ping among the frontier molecular orbitals. Due to their 
special structural motif, in K-phases the intradimer inter- 
action is larger than all interdimer interactions. There- 
fore, as already suggested by other authors @||, the 
bonding and antibonding dimer orbitals are a convenient 
basis for tight-binding calculations. 

The hopping integrals have been recently estimated for 
ET-Br by ab initio calculations . The intradimer hop- 
ping 6i=0.272 eV is more than two times larger than the 
interdimer integrals, p =0.13 eV and q = 0.04 eV along 
c + a, and 62 =0.085 eV along c direction. The dimer 
bonding and antibonding orbitals are split by an amount 
ss 2&i, whereas interdimer hoppings give residual inter- 
actions (p ± g)/4 and 62/2. By neglecting the mixing 
between bonding and antibonding orbitals, the original 
four site problem reduces to two independent two site 
problems. Since each dimer bears 3 electrons, the upper 
(conduction) band is half-filled. Focusing on antibonding 
orbitals, each dimeric site has 2 nearest neighbors along c 
direction (interaction 62/2), and 4 next-nearest neighbors 
along c ± a directions (interaction (p + q)/2). A funda- 
mental advantage of the dimer model is that the dimer 
basis fully exploits the local dimer symmetry, so that the 
resulting tight-binding model is more symmetric than 
the underlying lattice. The unit cell can therefore be re- 
duced to contain a single site. The analytical expression 
for the conduction band in the doubled Brillouin zone is: 

e(k) = 62 cos k z c + 2(p + q) cos cos 

In fig. 1 we report the energy dispersion curves for 
the conduction and valence bands, the density of states 
(DOS) for the conduction band, and the Fermi surface 
(FS), calculated in the dimer approximation, as com- 
pared with those obtained by solving the tight-binding 
problem with four interacting site orbitals. The main 
features of the complete four orbital model are well re- 
produced in the dimer approximation, namely, the band- 
width, chemical potential, FS topology and DOS, and 
also the location of the van Hove singularity (vHs). The 
bands in fig. 1 compare well with EHT bands ||, justi- 
fying a posteriori the neglect of inner orbitals. Recently 
published LDA bands for k-ET 2 Cu(NCS) 2 and ET-Br 
Q are also similar apart from an overall bandwidth nar- 
rowing. Moreover the topology of the FS agrees with 
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the experimental Shubnikov-de Haas (SdH) data on k- 
ET 2 Cu(NCS) 2 § and on k-ET 2 I 3 §. Therefore the 
dimer model offers a good description of the K-phase 
band structure. We adopt this model to investigate SC. 

Since in OSC the coherence length is short (of the order 
of few lattice spacings) [El, we choose a real space pair- 
ing potential analogous to that widely used for cuprate 
superconductors JlO| | : 

V = -J2^j a h a ]i a jl a ^- (O- 1 ) 



The sum is extended to nearest and next-nearest 
neighbors (with coefficents gNN and gNNN respectively), 
and also includes i — j (on-site) go interaction, with 
either positive (pairing) or negative (re puls ion) coeffi- 



cents. With the three parameters in Eq. |0.1| we explore 
all the possible symmetries of the gap function A(k) in 
the bidimcnsional point group pgg of the organic lay- 
ers of the crystal. In particular, the BCS gap func- 
tion for the singlet pairing can be written as a linear 
combination of functions with different symmetries [fl0|| : 
A(k) = Eti A ^( k X where m (k) = sin sin *f 
transforms as the B g irreducible representation, and 
7? 2 (k) = cos %r cos 773 (k) = 1, and 774 (k) = cosk z c 
as A g . Adopting the standard notation for cuprates 
771 corresponds to d xy pairing and 773 to sq isotropic pair- 
ing, while ?7 2 and 774 correspond to the extended or gener- 
alized s* pairing. The d xy symmetry in the K-phase OSC 
corresponds to d x i_ y i in the cuprates, due to the differ- 
ent orientation of the crystallographic axis with respect 
to the underlying lattice. We observe that, at variance 
with cuprates, d-wave cannot mix with s-wave compo- 
nents. 

Inserting the experimental T c = 11.6 K, relevant to 
ET-Br, into the gap equation, only two free pairing pa- 
rameters survive, that we choose as the ratios of go and 
gNN over gNNN- We solve the BCS gap equation on 
a 1888x1248 lattice in the go~gNN parameter space to 
obtain the phase diagram in fig. 2. For several points in 
the phase diagram the calculated gap has been checked 
to correspond to the absolute minimum of the BCS free 
energy. At T = T c the phase diagram is very simple: the 
gap equation is linear so that the mixing of pair functions 
with different symmetry is strictly forbidden. By lower- 
ing T below T c , the non-linearity of the BCS equation al- 
lows the mixing of s and d pairing, and the mixing region 
widens with decreasing T. However, the actual amount 
of mixing is always very small (< 10~ 3 ). In fact, T c is 
rather low, so that the gap equation stays quasi-linear, 
and the mixing of gap functions with different symmetry 
is small |yj . 

Having solved BCS gap equation, we calculate the 
macroscopic properties of the superconducting state us- 
ing standard approaches iQ. Since in the dispersion 
curves of the normal state the vHs occurs at energies 
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much higher than the superconducting gap (fig.l), the 
superconducting properties are largely dominated by the 
states on the FS. In particular, low energy macroscopic 
properties are strongly affected by the presence of nodes 
in the gap function at the FS. On the other hand the max- 
imum gap at the FS corresponds to a divergence in the su- 
perconducting DOS, and, as such, dominates frequency- 
dependent susceptibilities. In fig. 3 we report the mini- 
mum and maximum of |A(k)| calculated at T = 0. In 
the <i-wave region, the minimum gap at the FS vanishes, 
i.e. nodes are observed in the superconducting gap at the 
FS. In the s*-region pseudonodes (not necessarily imply- 
ing a change of sign of the gap at the FS), are observed. 
At small |<7aw| the s-wave gap is mainly isotropic and a 
sharp crossover is observed between d-wave and s-wave 
regions. Due to the negligible mixing of d— and s-pair 
functions, in the d-wave region the maximum gap is fixed 
by the experimental T c to 2.09 T c (see Fig. 3, bottom 
panel). In the s-region, due to the competition between 
isotropic and anisotropic s-wave components the maxi- 
mum gap at the FS shows a more complex structure. We 
observe that for very small jgjvjvl the gap is dominated 
by the isotropic component with the BCS value 1.76 T c . 

For the three representative points marked in the 
phase diagram in fig. 2, where the gap is dominated by 
the isotropic, anisotropic s-wave or d-wave contribu- 
tion, respectively, we calculate the superconducting elec- 
tronic specific heat, magnetic penetration depth, and 
tunneling bulk conductance. The low-energy behavior of 
these quantities is exponential-like in the isotropic case, 
whereas it is power-law in the anisotropic case, in agree- 
ment with standard results for conventional BCS Jl^ ] 
and rf-wav e JL3| SC. Experimental estimate of the spe- 
cific heat ]14| and tunneling |l5j are affected by large 
uncertainties, so that the comparison is not conclusive. 
On the other hand, experimental studies of the tempera- 
ture dependence of the magnetic penetration depth led to 
conflicting conclusions, suggesting either isotropic pair- 
ing or anisotropic gapless SC Jl^|. The superconducting 
phonon self-energy, presents a singularity at a frequency 
corresponding to twice the maximum of the gap at the 
FS |T7j. By lowering T below T c , phonons lying below 
(above) this frequency soften (harden). The crossover 
frequency strongly depends on the topology of the FS, so 
that the choice of the model for the electronic structure 
is now crucial. The calculated crossover of the phonon 
self-energy does in fact coincide with 2max|A(k)|. As 
shown in fig. 3, max|A(k)| reaches its minimum value 
for isotropic So pairing, corresponding to the BCS value 
2Ao = 3.53T C as 28 cm _1 . The salient feature of a recent 
low-energy Raman scattering study jl8| is that below T c 
all relevant bands harden with the greatest relative shift 
exhibited by the lowest observed phonon at 27.4 cm -1 . 
This value corresponds to a lower limit for the maximum 
of the gap at the FS. On the basis of our model, this 
result can only be interpreted assuming isotropic so pair- 
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ing. Similar results hold for k-(ET) 2 Cu(NCS) 2 , as can be 
inferred from recent neutron scattering experiment Jl9| . 

The presented picture for SC in re-phases is based on 
a single-particle description of the electronic structure. 
We have extracted the single-particle parameters from 
HF-SCF calculations 0: the comparison with EHT @ 
and LDA || bands, as well as with experimental data 
on k-(ET) 2 Cu(NCS) 2 [@ and on re-(ET) 2 I 3 § indicates 
that the resulting picture is extremely accurate in re- 
producing the topology of the FS. On the other hand, 
residual many-body effects are important in determin- 
ing susceptibilities or dynamical properties of OSC 
For instance, the calculated electronic specific heat is 
quasi-linear with T, but the slope, 7 ~ 7.5 mJ/mole ■ 
K 2 , is three times smaller than the experimental value, 
22 ± 3 mJ/mole-K 2 for ET-Br or 24 ± 3 mJ/mole-K 2 for 
k-(ET) 2 Cu(NCS) 2 0. Moreover, the effective masses 
inferred from SdH data via the Lifshits-Kosevich formula 
for k-ET 2 Cu(NCS) 2 § and k-ET 2 I 3 § are about three 
times larger than the single particle effective masses. It is 
important to stress that our picture for SC survives also 
in the presence of residual many-body interactions, in 
the hypothesis that these interactions do not modify the 
topology of the FS, but only lead to an overall renormal- 
ization of single- particle parameters. With the chosen 
pairing potential, this only implies an overall rescaling of 
pairing parameters: for instance, by rescaling our band 
parameters by a factor of 2, the same phase diagram as 
in fig. 2 is obtained, with pairing interactions reduced by 
a factor close to 2. 

On the other hand, the large variety of physical prop- 
erties characterizing K-phases suggests that residual in- 
teractions can also act to modify the FS topology. In par- 
ticular, ET-C1 is an antiferromagnetic insulator, whereas 
SdH oscillations are not observed in the (conducting) 
Br-analog at ambient pressure. High pressure SdH os- 
cillations in ET-Br point to a distortion of the FS [po| . 
In this connection, the proposed dimer model is simple 
enough to offer clues about the fine tuning of the physical 
properties of K-phases as induced by variations of the FS 
topology. By looking at the band-structure in fig.l, one 
realizes that an extended vHs is present along the MZ 
direction, at energy — 6 2 with respect to the Fermi level. 
Even a small reduction of 6 2 would imply a large increase 
of the effective mass (due to the large effect of the vHs), 
and a reduction of the a orbit, in agreement with SdH 
measurement in ET-Br [ pof . In the extreme limit of van- 
ishing fe 2 the vHS would lie at the FS, and the FS itself 
would be a rectangle with perfect nesting. This extreme 
picture seems expecially attractive to understand AFM 
in ET-C1 ||] as well as the competition between SC and 
AFM in ET-C1 and in ET-Br. 6 2 is a small parameter 
and therefore it is subject to rather large relative vari- 
ations in the different structures. In Et-Cl and ET-Br 
salts the CI and Br atoms are located very near to the 6 2 
dimer 0, possibly affecting the corresponding transfer 
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integral. Ab initio calculations are in progress to test the 
effect of the counterions and of residual electron-electron 
interactions on the 62 hopping. 

A different source of renormalization of the hopping 
parameters involves electron-phonon interactions. As it 
is evident from fig.l, there are large regions in the k space 
where the effective electronic masses are very large, even 
larger than the phonon masses. In other words, since 
the ab initio &2 value corresponds to 685 cm -1 , so that 
most of the intramolecular phonons |^l|,^2| arc in the an- 
tiadiabatic limit with respect to 62 hopping. Along this 
direction an antiadiabatic polaron narrowing j2^] can be- 
come effective, reducing 62 hopping. From the available 
coupling constants and phonon frequencies ]2l]]2^] , the 
polaron narrowing factor is estimated about 3. 

In summary, we have shown that an extremely simple 
model - a tight-binding lattice with one site per cell - 
offers an accurate description of the band structure of 
K-phase ET salts. Using this model we were able to get 
a solution of the full, non-linear BCS gap equation for 
OSC. At variance with cuprates, our simple one-band 
model can be derived from first-principle calculations, 
and therefore allows us to get a significant comparison 
with experiment. Our work supports isotropic sq SC. 
Finally, the possible role of counterions as well as of 
phonon-mediated interactions is invoked to justify the 
large variability of physical properties of K-phase ET 
salts. 
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Figure captions 



Figure 1. The tight-binding highest occupied bands (left 
panel), DOS (right panel), and FS (inset) calculated for 
K-(ET) 2 Cu[N(CN) 2 ]Br in the dimer approximation (solid 
line), and for the full four site problem (dotted line). 

Figure 2. The three-dimensional phase diagram for SC. 
Dashed lines are the isotherms for temperatures below 
T c . Points lying between two isotherms correspond to 
mixed s- and d-wave states. The circle, diamond, and 
star indicate representative points for isotropic (gNNN — 
4.96 • 10 -2 eV), anisotropic s-wave (gNNN — 4.76 • 10~ 2 
eV), and d xy -wa,ve (gNNN — 7.24- 1CT 2 cV), respectively. 

Figure 3. The minimum (top panel) and maximum (bot- 
tom) of |A(k)| on the Fermi Surface at T = 0, as a 
function of the relative pairing parameters strengths (see 
text). 
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